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METHOD OF PROCES SING MARINE MAGNETIC GRADIENT DATA awn 
EXPLORATION METHODS USING THAT DATA 

Field of the Invention 

This invention relates to a method of processing marine 
magnetic gradient data and also to exploration methods 
using that data. 

Background Art 

Exploration of the earth to locate oil, minerals or other 
valuable deposits, uses seismic techniques in order to 
provide information about the subsurface structure of the 
earth so that an assessment can be made as to whether a 
particular deposit may or may not be present. in 
particular, seismic techniques provide an indication of 
various subsurface structures including porous sandstone 
or fractured carbonates which may contain oil or other 
valuable deposits. 



Whilst this technique does provide an indication of 
subsurface structures, it is possible a particular 
structure which may appear to be of interest, in fact 
represents a significant hazard to drilling. For example, 
it is possible that subsurface volcanos exhibit very 
similar characteristics to anticlines which may contain 
oil. if a drilling operation is initiated and it is found 
that the anticline is in fact a volcano, then a 
considerable financial loss results because of the cost of 
establishing the drilling platform. This is particularly 
the case in relation to marine exploration because of the 
significantly higher drilling cost. 

in order to provide more information concerning the nature 
of subsurface structures, magnetic data of the exploration 
3 5 area is obtained. 
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A magnetic gradient survey enables magnetic data relating 
to the survey region to be obtained and this magnetic data 
can be used to provide information relating to the nature 
of subsurface structures. if the magnetic survey is 
overlayed with the seismic survey, then structures which 
appear to be of interest from the point of view of the 
seismic survey, can be further considered in the light of 
the magnetic data so that a clearer indication can be 
formed as to whether the structure is a structure which 
may contain a deposit of interest such as an oil deposit, 
or whether the structure exhibits magnetic phenomena, such 
would be the case with a volcano. Thus, the location of 
drilling platforms can be decided with more precision to 
avoid subsurface structures which could be hazardous from 
15 the point of view of a drilling operation. However, 

conventional techniques for processing magnetic gradient 
data contain considerable distortion and anomalies due to 
unwanted magnetic effects including the ship bias referred 
to above. 

20 

The conventional method of obtaining that data in relation 
to onshore exploration is to tow magnetometers behind an 
aircraft to obtain magnetic data in relation to the survey 
area. In the case of marine exploration, the 
25 magnetometers are towed behind a ship. 



The usefulness of marine magnetic data which is obtained 
in this manner has been limited by the quality of the data 
which is obtained. One of the major problems with marine 
30 magnetic data acquisition is interference or so-called 

ship bias which is created by the magnetic field induced 
by the ship which tows the magnetometers. 



In conventional marine magnetic data acquisition, two 
magnetic field sensors, which are generally referred to 
fish, are towed behind a ship. The magnetometers are 
connected to a tow line and the magnetometer which is 
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closest to the ship is towed at a distance of some 300-600 
metres behind the ship to avoid the effect of ship induced 
magnetic field. Furthermore, the magnetometers are 
separated by a distance of in excess of 100 metres. The 
reason for the length of the tow line and the separation 
of the magnetometers is to reduce ship bias and therefore 
provide data which is relatively free of that bias. 

However, because the towing line is so long and the 
distance of separation so great, the degree of drift of 
the magnetometers in the ocean as the magnetometers are 
towed behind the ship is considerable. Furthermore, the 
processing techniques used to acquire magnetic gradient 
data assume that after a given time period, the trailing 
magnetometer will be towed to a position which coincides 
with a previous position of the magnetometer which is 
closest to the ship. This assumption is made during the 
processing of the magnetic data. However, because the tow 
line is extremely long and the distance of separation 
between the magnetometers is considerable, drift of the 
magnetometers due to sea currents and the like means it is 
unlikely that the trailing magnetometer will actually 
occupy the same position as the first magnetometer after 
that given time period. 

Furthermore, in the acquisition of the magnetic data, the 
ship is required to travel along predetermined survey 
lines and because of the drifting of the sensors, the 
sensors do not actually travel along the survey line. in 
the southern hemisphere, if the magnetometers drift to the 
north of the survey line, the recorded gradient data 
between the two magnetometers will have larger distortion 
than if drifting towards the south of the survey line. 

3 5 Summary of the Invention 

The object of the invention is to provide an improved 
processing method which accounts for ship-induced bias and 
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which therefore enables the magnetometers to be towed by a 
tow line which is much shorter than the conventional tow 
line and also be located closer together to reduce the 
amount of drift and also increase the likelihood that the 
5 trailing magnetometer will occupy the same position which 
was previously occupied by the leading magnetometer after 
a given time period. 

The invention may be said to reside in a method of 
processing marine magnetic data obtained by towing first 
and second spaced apart sensors behind a ship," the method 
including: 

obtaining raw magnetic gradient data from the 

sensors; 

determining the trend of the gradient of the ship 
bias detected by the sensors; 

subtracting the trend from the raw magnetic 
gradient data to obtain corrected gradient data; and 

processing the corrected gradient data to provide 
a data output. 

Because the method determines the trend of the gradient of 
the ship bias, a more accurate estimation of the actual 
gradient bias introduced by the ship, having regard to any 
25 drifting of the sensors, is taken into account. By 

subtracting the trend of the gradient bias from the raw 
gradient data, the bias introduced by the ship can be 
removed from the data so that more accurate data is 
obtained. Because the present method eliminates the 
30 effect of the ship bias, the sensors can be towed much 
closer to the ship and at much closer separation 
distances, thereby reducing the amount of drift of the 
sensors and therefore resulting in the sensors being more 
likely to follow the actual survey line upon which the 
35 ship travels and also for the rear sensor to occupy the 
same position at a later time which was previously 
occupied by the front sensor. Thus, the quality of the 
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output data is greatly improved and therefore provides 
data which gives a much more accurate indication of the 
subsurface magnetic characteristics of the survey area, 
which can be then used in combination with seismic data to 
5 make an assessment as to the viability of further 
exploration or a drilling, program. 

Preferably the method includes the step of determining an 
estimate of the gradient of the ship bias from the raw 
10 magnetic gradient data obtained by the sensors, and 

determining the trend of the gradient of the ship bias 
from that estimate of the gradient of the ship bias. 

Preferably the raw magnetic gradient data is calculated as 
15 the measured magnetic signal at the leading sensor less 
the measured magnetic signal at the trailing sensor, 
divided by the distance between the sensors. 

The estimate of the gradient of the ship bias may in fact 
20 simply be derived from the raw magnetic gradient data 

which is obtained from the sensors and which may include 
data relating to the gradient of the ship bias as well as 
data relating to other magnetic effects. 

25 Preferably the trend of the gradient of the ship bias is 
determined by applying the gradient of the ship bias data 
to a B-Spline filter with smoothness 0.8. 

Preferably the processing comprises the application of the 
30 corrected gradient data over sample intervals, integration 
of the corrected gradient data into total magnetic 
intensity data, and application of a low pass filter to 
the integrated total magnetic intensity data to provide 
the data output. 

35 

Preferably the total magnetic intensity is obtained by 
integrating the corrected gradient data in accordance with 
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the following equation: 

TMI=£ [G c Ax f (0] + M 5WI 
/=! 

wherein G c is the corrected gradient data 
5 obtained after subtraction of the trend of the gradient of 
the ship bias from the raw gradient data, and 

M s tat is the total magnetic intensity value at a 
point representing the start of the survey or at a 
location where a survey line and a tie line cross each 
10 other. 

The values, of TMI so calculated are then smoothed to 
remove residual noise using a low pass filter • 

15 Preferably the gradient of the ship bias is calculated 
according to the following equation: 

Gbias = M / W " M - W , Wt,)-D{t 2 )] + {M b {t,)-M h {t 2 K 
A/ A/ 

2 0 where : 

Mj(x) = M e (x) + D(t x ) + M b (t x ) 
Mfc) = M e (jc) + D{h) + M b (t 2 ) 

25 

where M f is the magnetic field as measured by the 
front sensor and which consists of the environmental field 
Me(x), diurnal variation D(ti), ship bias M b (t x ) caused by 
ship-induced field, sensor drifting and sensor 

3 0 perturbation, all at time t x and along track distance x # 

and at some later time t 2 , Al is the distance between the 
sensors, D<t 2 ) is the diurnal variations sensed by the rear 
sensor, and 

M b (t 2 ) is the ship bias at time t 2 detected by the 
35 rear sensor M r . 
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In one embodiment of the invention the first ana second 
sensors which are towed behind the ship are included in a 
group of three or more towed sensors. 

5 In this embodiment, preferably the number of sensors towed 
behind the ship comprises three sensors. 

According to this embodiment data from any two of the 
sensors is used to provide the raw magnetic gradient data. 

10 

However, the embodiment may obtain raw gradient data from 
all three sensors and determine the trend of the gradient 
of the ship bias detected by all three sensors. 

15 The invention may also be said to reside in a method of 
obtaining gradient data for an exploration program, the 
method including: 

towing first and second sensors behind a ship 
along predetermined survey lines; 

20 obtaining raw magnetic gradient data from the 

sensors ; 

determining the trend of the gradient of the ship 
bias detected by the sensors; 

subtracting the trend from the raw magnetic 
25 gradient data to obtain corrected gradient data; and 

processing the corrected gradient data to provide 
a data output. 



Preferably the method includes the step of determining 
30 estimate of the gradient of the ship bias from the raw 
magnetic gradient data obtained by the sensors, and 
determining the trend of the gradient of the ship bias 
from that estimate of the gradient of the ship bias. 

35 Preferably the raw magnetic gradient data is calculated 
the measured magnetic signal at the leading sensor less 
the measured magnetic signal at the trailing sensor, 
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divided by the distance between the sensors. 

The estimate of the gradient of the ship bias may in fact 
simply be derived from the raw magnetic gradient data 
5 which is obtained from the sensors and which may include 
data relating to the gradient of the ship bias as well as 
data relating to other magnetic effects. 

Preferably the trend of the gradient of the ship bias is 
10 determined by applying the gradient of the ship bias data 
to a B-Spline filter with smoothness 0.8. ■ 

Preferably, the further processing comprises the 
application of the corrected gradient data over sample 
15 intervals, integration of the corrected gradient data into 
total magnetic intensity data, and application of a low 
passed filter to the integrated total magnetic intensity 
data to provide the data output. 

20 Preferably the total magnetic intensity is obtained by 

integrating the corrected gradient data in accordance with 
the following equation: 

TMI=£ [GcAx.it)} +M sla , 
i=i 

25 wherein G c is the corrected gradient data 

obtained after subtraction of the ship bias trend from the 

raw gradient data, and 

M„tat is the total magnetic intensity value at a 

point representing the start of the survey or at a 
30 location where a survey line and a tie line cross each 

other . 

The values of TMI so calculated are then smoothed to 
remove residual noise using a low pass filter. 

35 

Preferably the gradient of the ship bias is calculated 
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according to the following equation: 

Gbias = M 'M- M **K [£>(0-Z>(r 2 )] + [M,(/,)-A/ fc (r ? )] 



where : 



Mj(x) = M e (x) + D(t { ) + M b {U) 
M£x) = M e (x) + D(t 2 ) + M b (t 2 ) 



where M f is the magnetic field as measured by the 
front sensor and which consists of the environmental field 
Me(x), diurnal variation D(tx), ship bias M b (^) caused by 
ship-induced field, sensor drifting and sensor 
15 perturbation, all at time tl and along track distance x, 
and at some later time t 2 , Al is the distance between the 
sensors, D(t 2 > is the diurnal variations sensed by the rear 
sensor, and 

M b (t 2 ) is the ship bias at time t 2 detected by the 
2 0 rear sensor M x . 

In one embodiment of the invention the first and second 
sensors which are towed behind the ship are included in a 
group of three or more towed sensors. 



In this embodiment, preferably the number of sensors towed 
behind the ship comprises three sensors. 

According to this embodiment data from any two of the 
sensors is used to provide the raw magnetic gradient data. 

However, the embodiment may obtain raw gradient data from 
all three sensors and determine the trend of the gradient 
of the ship bias detected by all three sensors. 

The invention may also be said to reside in a method of 
subsurface exploration to determine viability of drilling 
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in a marine environment by considering magnetic data 
relating to the environment, and wherein the magnetic data 
has been obtained by: 

obtaining raw magnetic gradient data from the 

sensors ; 

determining the trend of the gradient of the ship 
bias detected by the sensors; 

subtracting the trend from the raw magnetic 
gradient data to obtain corrected gradient data; and 

processing the corrected gradient data to provide 
a data output. 

Preferably, obtaining the magnetic data includes the step 
of determining the gradient of the ship bias from data 
obtained by the sensors, and determining the trend of the 
gradient of the ship bias from the gradient of the ship 
bias. 

Preferably the raw magnetic gradient data is calculated as 
the measured magnetic signal at the leading sensor less 
the measured magnetic signal at the trailing sensor, 
divided by the distance between the sensors. 

The estimate of the gradient of the ship bias may in fact 
25 simply be derived from the raw magnetic gradient data 

which is obtained from the sensors and which may include 
data relating to the gradient of the ship bias as well as 
data relating to other magnetic effects 
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Preferably the trend of the gradient of the ship bias is 
determined by applying the gradient of the ship bias data 
to a B-Spline filter with smoothness 0.8. 

Preferably the processing comprises the application of the 
corrected gradient data over sample intervals, integration 
of the corrected gradient data into total magnetic 
intensity data, and application of a low passed filter to 
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the integrated total magnetic intensity data to provide 
the data output. 

Preferably the total magnetic intensity is obtained by 
integrating the corrected gradient data in accordance with 
the following equation: 

TMI=£ [Gc&xfr)] + M sla , 
i=i 

wherein G c is the corrected gradient data 
obtained after subtraction of the ship bias trend from the 
raw gradient data, and 

M stat is the total magnetic intensity value at a 
point representing the start of the survey or at a 
location where a survey line and a tie line cross each 
15 other. 

The values of TMI so calculated are then smoothed to 
remove residual noise using a low pass filter. 

20 Preferably the gradient of the ship bias is calculated 
according to the following equation: 

Gbias = M/M-JM*), [DO, )-£>(/,)] + [M h 



2 5 where : 



30 



Mf(x) = M e (x) + £(/,) + M b (h) 
M£x) = M e (x) + D(t 2 ) + M b (t 2 ) 



where M £ is the magnetic field as measured by the 
front sensor and which consists of the environmental field 
Me(x), diurnal variation D(ti), ship bias M b (t!) caused by 
ship-induced field, sensor drifting and sensor 
3 5 perturbation, all at time t x and along track distance x, 
and at some later time t 2 , Al is the distance between the 
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sensors, D(t a ) is the diurnal variations sensed by the rear 
sensor, and 

M b (t 2 ) is the ship bias at time t a detected by the 
rear sensor M r . 

In one embodiment of the invention the first and second 
sensors which are towed behind the ship are included in a 
group of three or more towed sensors. 

In this embodiment, preferably the number of sensors towed 
behind the ship comprises three sensors. 

According to this embodiment data from any two of the 
sensors is used to provide the raw magnetic gradient data. 

However, the embodiment may obtain raw gradient data from 
all three sensors and determine the trend of the gradient 
of the ship bias detected by all three sensors. 



The invention may still further be said to reside in a 
method of drilling for a deposit in a marine environment, 
including: 

determining the location of drilling from data 
which has been obtained and which indicates the possible 
25 existence of the deposit; and 

which location is also determined by magnetic 
data which has been obtained by towing magnetic sensors 
behind a ship, the magnetic data being processed by: 

obtaining raw magnetic gradient data from the 
3 0 sensors ; 

determining the trend of the gradient of the ship 
bias detected by the sensors; 
subtracting the trend from the raw magnetic 
gradient data to obtain corrected gradient data; 
35 and 

processing the corrected gradient data to provide 
a data output. 
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Preferably the method includes the step of determining an 
estimate of the gradient of the ship bias from the raw 
magnetic gradient data obtained by the sensors, and 
5 determining the trend of the gradient of the ship bias 
from that estimate of the gradient of the ship bias. 

Preferably the raw magnetic gradient data is calculated as 
the measured magnetic signal at the leading sensor less 
the measured magnetic signal at the trailing sensor, 
divided by the distance between the sensors. 

The estimate of the gradient of the ship bias may in fact 
simply be derived from the raw magnetic gradient data 
which is obtained from the sensors and which may include 
data relating to the gradient of the ship bias as well as 
data relating to other magnetic effects. 

Preferably the trend of the gradient of the ship bias is 
determined by applying the gradient of the ship bias data 
to a B-Spline filter with smoothness 0.8. 

Preferably the processing of the corrected gradient data 
comprises the application of the corrected gradient data 
over sample intervals, integration of the corrected 
gradient data into total magnetic intensity data, and 
application of a low passed filter to the integrated total 
magnetic intensity data to provide the data output. 

Preferably the total magnetic intensity is obtained by 
integrating the corrected gradient data in accordance with 
the following equation: 

TMI=£ [C c Ax,(/))+A^ 



*stm 



35 wherein G c is the corrected gradient data 

obtained after subtraction of the ship bias trend from the 
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raw gradient data, and 

M 8tat is the total magnetic intensity value at a 
point representing the start of the survey or at a 
location where a survey line and a tie line cross each 
5 other. 

The values of TMI so calculated are then smoothed to 
remove residual noise using a low pass filter. 

10 Preferably the gradient of the ship bias is calculated 
according to the following equation: 

= "/(*)-*,(*) = [DiV-D^ + lMM-MM) 
A/ Al 

1 5 where : 

Mj(x) = M e (x) + £>(*,) + M b (t x ) 
M£x) = M e (x) + D(t 2 ) + M b (t 2 ) 

where M £ is the magnetic field as measured by the 
front sensor and which consists of the environmental field 
M*(x), diurnal variation D(t x ) , ship bias M b ( tl ) caused by 
ship-induced field, sensor drifting and sensor 
perturbation, all at time tl and along track distance x, 
and at some later time t 2 , Al is the distance between the 
sensors, D(t 2 ) is the diurnal variations sensed by the rear 
sensor, and 

M b (t 2 ) is the ship bias at time t 2 detected by the 
30 rear sensor M r . 

In one embodiment of the invention the first and second 
sensors which are towed behind the ship are included in a 
group of three or more towed sensors. 
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In this embodiment, preferably the number of sensors towed 
behind the ship comprises three sensors. 
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According to this embodiment data from any two of the 
sensors is used to provide the raw magnetic gradient data. 

However, the embodiment may obtain raw gradient data from 
all three sensors and determine the trend of the gradient 
of the ship bias detected by all three sensors. 

Brief Description of the Drawings 

A preferred embodiment of the invention will be described, 
by way of example, with reference to the accompanying 
drawings in which: 

Figure 1 is a view showing a marine magnetic data 

survey; 

Figure 2 is a flow chart showing the preferred 
method of processing the marine magnetic gradient data 
according to the preferred embodiment of the invention; 

Figure 2A is a graph showing the function of a 
low pass filter used in the preferred embodiment; 

Figure 3A and Figure 3B are graphs showing data 
obtained according to the survey and as processed 
according to the preferred embodiment of the invention; 

Figure 4 is a representation of a survey using 
data created according to a prior art technique; 

Figure 5 is a diagram similar to Figure 4, but 
using data obtained according to the preferred embodiment 
of the present invention; and 

Figure 6 is a view of a further embodiment of the 

invention. 



Description of the Preferred Embodiment 

With reference to Figure 1, a marine survey to acquire 
magnetic data is illustrated. A ship 10 tows first 
magnetometer Ml and second magnetometer M2 by a tow line 
35 12. The ship travels along survey line 14 and the line 

16, which are dotted, represent the ship-induced magnetic 
field. 
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In Figure 1, it is assumed the ship is in the southern 
hemisphere and is heading north-northwest, and the sensors 
are drifting towards the northeast. In this case, the 
recorded gradient data between sensors Ml and M2 will have 
a larger distortion than if the sensors were drifting 
towards the southwest. 

Ideal conditions would assume that the sensors Ml and M2 
are towed directly behind the ship on the survey line and 
that therefore, after a given time period, the- sensor M2 
would occupy the same position as was previously occupied 
by the sensor Ml. However, due to drifting of the sensors 
because of ocean currents and the like, this ideal 
15 situation in general cannot be achieved. Nevertheless, 
the processing methods used to process magnetic data do 
assume that the sensor M2 will occupy the same position as 
the sensor Ml, previously occupied, but at a later time. 
The processing method according to the preferred 
embodiment enables ship bias to be accounted for much more 
accurately than any prior art techniques and therefore 
enables the sensors to be towed closer to the ship and to 
be separated by a smaller distance. This therefore 
provides much greater accuracy in the measurements because 
the amount of drift is not as great due to the shorter tow 
line and also because the separation between the sensors 
is smaller, the likelihood of the rear sensor occupying 
the same position as the front sensor at a later period in 
time is greatly increased. 
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At any time t, and along-track distance x, the front sensor 
measures magnetic field - Mf consisting of the environment 
field Me(x), the diurnal variation U(t,), ship bias t^{t x ) 
caused by ship induced field, sensor drifting and sensor 
3 5 perturbation . 



Mj(x) = M e (.x) + D(t,) + M b (ti) 
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At some later time t 2 , the rear sensor makes a measurement 
at the same spatial point: 

Mr(x) = M e (x) + D(t 2 ) + M b (t 2 ) 

The difference between the two measurements at the same 
location is: 

M{x) - MXx) = [D(t0 - D(t 2 )] + [M b ( tl ) - M b {t 2 )) 
Ship bias gradient (Step 2, Figure 2) can be expressed as: 

Gbias = U>(/,)-Z)(/,)] + rM,(f,)-M^)j 
A/ Xl : 



It should be mentioned that the calculated ship bias 
gradient is still affected by diurnal variation. in 
practice we notice that the sea current drifting is a kind 
of long wave-length variation so ship bias should be 
coincident with such current drifting. Step 2 in Figure 
2, which is the calculation of the gradient of the ship 
bias, is actually an estimate of the gradient of the ship 
bias which is obtained from the raw data collected by the 
sensors Ml and M2 . The data collected by the sensors will 
25 include raw magnetic gradient data which contains many 

signals including environmental magnetic signals, diurnal 
signals and the ship bias, as well as instrument bias and 
drift. The trend of the gradient of the ship bias is 
obtained from this estimate by applying the estimate of 
the gradient of the ship bias to a B-Spline filter, as 
will be described in more detail hereinafter. Therefore 
the trend of the gradient of the ship bias (Gtrend) (Step 
3, Figure 2) is used for the correction along a track 
line. 



The gradient correction of the ship bias (<? f )can be 
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expressed as follows: 

G C = G- G trt nd (Step 4, Figure 2) 
5 in this equation, G is the raw magnetic gradient data. 

As said before such drifting can cause large bias effects. 
Consequently, when the sensor drifting and perturbation 
take place simultaneously the bias zig-zag will be around 
10 its mean values. The variable mean value of the bias 

along a survey line is regarded as bias trend Isince the 
high frequencies of sensor perturbation only create random 
noise around the bias trend and after integration its 
effect can be eliminated from raw gradient data. 
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After removal of ship bias effects TMI (Total Magnetic 
intensity) can be calculated (Step 6) through integration 
of magnetic gradient data: 



20 TMI= ■•=■ [GcAxi(t)]+Mstat 

Here, Ax^t) is the sampling distance along the survey 
line. W gtat is the TMI value at a point of survey start or 
at the location where the survey line and a tie line cross 
25 each other. 

The TMI value calculated in Step 6 is then smoothed using 
a low pass filter (Step 7) to remove any features which 
have a greater rate of change of TMI with distance than is 
expected in the particular survey area. An example of the 
action of this filter is provided in Figure 2A in which 
trace 50 is the smooth TMI curve and trace 51 is the TMI 
data prior to smoothing. 

35 The data output obtained at Step 7 may include line 

levelling and data griding to provide a final data output. 
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Figures 3A and 3B represent an actual example of the 
preferred embodiment of the invention which relates to raw 
data obtained from a known region. 

In Figure 3A, the trace 20 represents the total ship bias. 
The line 21 represents the bias trend and it can be seen 
from the left hand side of the trace 20 to the right hand 
side of the trace 20, the bias trend 21 changes non- 
uniform^ about the value of -0.08 on the graph of Figure 
3A. A particular value for the bias trend can therefore 
be obtained for various intervals and this value can be 
subtracted from the raw data in order to provide the 
corrected data. The trend of the gradient of the ship 
bias is preferably determined by calculating the gradient 
of the bias data 20 and then applying that gradient data 
to a B-Spline filter with smoothness 0.8 in order to 
provide the representation of the ship bias 21. 

In the preferred embodiment of the invention, the gradient 
of the bias data 20 is calculated simply by using the raw 
magnetic gradient data obtained from the sensors because 
when that data is applied to the B-Spline filter, all that 
remains is the gradient of the ship bias component. The 
ship bias 21 can then be subtracted from the calculated 
raw gradient. 



The trend of the gradient of the ship bias is a non-linear 
function which is represented by the line 21 in Figure 3A. 
The line is a measure of how the bias changes over time as 
the sensors are towed behind the ship. As is clearly seen 
from line 21, the trend is not constant or merely an 
average, but rather represents the fluctuation of the bias 
gradient and at some times, is above the value of 0.08 and 
35 sometimes below that value. As noted above, the trend is 
determined by applying the gradient of the ship bias to a 
B-Spline filter with smoothness 0.8. However, in other 



WO 03/052460 PCT/AU02/01669 

- 20 - 

embodiments, the filter could have different smoothness 
depending on the region from which data is collected and 
the nature of the data which is collected. in general, 
the function of the filter is to smooth the trace 20 so 

5 that some meaningful value of the bias trend at particular 
time periods can be obtained. Thus, the filter is 
effectively determining the peaks and troughs of the trace 
20 and fitting a curved line between those peaks and 
troughs, which gives a measure of how the gradient of the 

.0 ship bias is changing with time as the sensors, are towed 
behind the ship. 



In the example of data obtained from the known region, 
dashed line 30 represents total magnetic intensity data 
which is provided according to the prior art processing 
technique. The line 32 represents the integrated total 
magnetic intensity data from the bias corrected gradient 
data according to the preferred embodiment of the 
invention. The line 34 represents the diurnal variation 
at a station some 500 km away from the survey area and 
line 36 is the observed field data including the diurnal 
effect. 



It can be seen from the conventional processing technique, 
which gives line 30, false anomalies are provided which do 
not show up on the line 32, which is produced according to 
the present invention. 



In general, the results of the magnetic survey are 
• 3.0 produced in a colour chart. Black and white 

representations of a colour chart are shown as Figure 4 in 
relation to conventional processing of data obtained from 
a known region and Figure 5 by a processing method 
according to the preferred embodiment of this invention. 



The characteristics of the Igneous basement of the known 
region are well known and it can be seen that data 
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processed according to the present invention produces a 
much clearer indication of actual magnetic structures than 
prior art technique which include considerable 
interference and data which may mislead analysers into 
5 considering that magnetic structures exist which are not 
actually there, or that magnetic structures do not exist 
when a magnetic structure actually is present. 
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Thus, the preferred embodiment of the present invention 
provides data which more accurately reflects the 
likelihood of magnetic structures which can then be used 
as a basis for determining the viability of a drilling 
operation in an exploration or recovery program. 

Figure 6 shows a second embodiment of the invention in 
which like reference numerals indicate like components to 
those described with reference to Figure 1. in this 
embodiment, three magnetometers Mi, M 2 and M 3 are towed by 
the ship 10 on tow line 12. This embodiment therefore 
builds some redundancy into the system in that, should one 
of the sensors fails, two sensors are still available in 
order to provide the required gradient data, and therefore 
there is considerably less likelihood that a survey will 
be wasted should a ship complete the survey only to find 
that one of the magnetometers has not operated properly. 
In this embodiment, preferably the distance between the 
magnetometers Mi and M 2 , and the distance between the 
magnetometers M 2 and M 3 is about 15 metres. The towing 
distance between the ship 10 and the first magnetometer Mi 
is preferably about 150 metres or less. 

This embodiment provides the further advantage that any 
group of two sensors can be used to provide the gradient 
data, and therefore the gradient data could be provided by 
the magnetometers M x and M 2 , the magnetometers M 2 and M 3 , or 
the magnetometers M x and M 3 . This embodiment provides the 
further advantage that all three magnetometers could be 
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used to provide the data to enable the calculation of the 
raw gradient and the gradient of the ship bias. The use 
of three magnetometers may improve accuracy because of the 
increased amount of data which is collected. 

5 

Obviously, more than three magnetometers could be used if 
desired, thereby increasing the amount of data collected 
and reducing the likelihood that a survey will be wasted 
because at least two of the magnetometers are not operated 
10 properly. 



Since modifications within the spirit and scope of the 
invention may readily be effected by persons skilled 
within the art, it is to be understood that this invention 
15 is not limited to the particular embodiment described by 
way of example hereinabove. 



